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We propose an approach for super-resolution optical lithography which is based on the inverse 
of magnetic resonance imaging (MRI). The technique uses atomic coherence in an ensemble of spin 
systems whose hnal state population can be optically detected. In principle, our method is capable 
of producing arbitrary one and two dimensional high-resolution patterns with high contrast. 
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I. INTRODUCTION 

Optical lithography is widely used to print circuit im¬ 
ages onto substrates because of its inherent parallelism 
and the fact that its non-ionizing illumination avoids sub¬ 
strate damage. However, the resolution of traditional 
optical lithoCTaphy is restricted by the Rayleigh diffrac¬ 
tion limit [il-lj. Due to this limit, we have to use shorter 
wavelengths such as extreme ultraviolet or X-ray to print 
smaller patterns i i- However, insulating materials 
such as silicon dioxide, when exposed to photons with 
energy greater than the band gap, release free electrons 
which subsequent can cause adverse charging which can 
lead to damage Q- It is therefore interesting and useful 
to invent an optical lithography scheme that can over¬ 
come the diffraction limit. 

In the past two decades, a number of methods have 
been proposed including m ultip hoton absor ption l8| -llfll|. 
quantum entanglement [lll - ll3l |. dopperlon |14l - ll6l| . and 
dark states El- However, these methods either require 
higher order multiphoton absorption, higher order quan¬ 
tum entangled states or multi atomic levels to write pat¬ 
terns with higher order resolution. Thus improving res¬ 
olution more than ten times the optical wavelength with 
these methods is difficult in practice. In 2010, we pro¬ 
posed a novel scheme to overcome the diffraction limit in 
optical lithography by inducing Rabi oscillations between 
two energy levels of the medium [Tsl - f^ . The extension 
from lower resolution to higher resolution is straightfor¬ 
ward where we just need to increase the power of the exci¬ 
tation laser or the interaction time. However, an effective 
way to print arbitrary one dimensional(lD) and two di¬ 
mensional (2D) patterns with low background is difficult 
in this method. Donut beam lithography .which is 

the inverse of the grou nd state depletion super-resolution 
(GSD) technique [^, can write arbitrary patterns 
with resolution beyond the diffraction limit, but it re¬ 


* prhemmer@ece.tamu.edu 

tzeyangliao@physics.tamu.edu 

tzubairy@physics.tamu.edu 


quires point-by-point scanning instead of full frame ex¬ 
posure and is very inefficient since most of the laser 
light is actually used to suppress resist exposure. Using 
frequency-encoded light masks in neutral atom lithogra¬ 
phy it is possible to print arbitrary sub-diffraction-limited 
features [^ . However, in atom lithography very high 
atom collimation is required for high resolution which is 
difficult to achieve with sufficient throughput to be useful 
for most lithography applications.. 

In this paper, we propose an approach based on the 
MRI super-resolution technique [^. Our method can 
be conveniently called magnetic resonance lithography 
(MRL) [i^. This method can be highly parallel, where 
the parallelism results from the fact that coherent excita¬ 
tion of two-level spin-flip transitions is employed. In this 
method, the extension from lower resolution to higher 
resolution is also straightforward where we just need to 
apply larger magnetic gradient and narrower bandwidth 
of microwave. Furthermore, it is also possible to print 
arbitrary ID and 2D patterns with low background. 

This paper is organized as follows. In Sec. H, we dis¬ 
cuss the basic physics of MRL and propose the pulse 
sequence to print arbitrary ID and 2D patterns. In Sec. 
HI, we numerically demonstrate how to print an arbitrary 
one dimensional and two dimensional patterns into the 
spin states and the photoresist. In Sec. IV, we discuss 
some possible materials for MRL. Finally, we summarize 
the results. 


II. BASIC PRINCIPLE AND PROPOSED 
EXPERIMENT 

MRL requires a paramagnetic system with high spin 
transition quality factors and high contrast optical read¬ 
out of the spin state. In addition, to write an arbitrary 
2D patterns the system should have three magnetic sub- 
levels. The proposed schematic setup is shown in Fig. 
1 (a) where the paramagnetic system is imbeded near the 
surface of the substrate. We can create the magnetic 
gradient fields using conductive metal wires fabricated 
on the substrate surface by means of a conventional pho¬ 
tolithography process. A static magnetic field Bq is ap- 
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FIG. 1: (Color online) (a) Schematic setup for magnetic reso¬ 
nance lithography. Bq is static magnetic field applied to split 
the magnetic sublevles m = ±1 while B(u) is microwave ap¬ 
plied to couple the transition m = 0 to m = ±1; (b) Energy 
levels for MRL; (c) The pulse sequence to write an arbitrary 
2D binary pattern, (d) An example showing how to write a 
letter “T” without crosstalk between lines. 


plied from the perpendicular direction to uniformly split 
the magnetic sublevels m = 0 and m = ±1 shown in 
Fig. 1(b). A microwave pulse with special spectral com¬ 
ponents is also applied from the perpendicular direction 
to excite the spin transitions m = 0v^m = ±l. Af¬ 
ter preparing the required pattern in the spin states, we 
may use a laser to excite the spins in state m = 0 into 
a fluorescing excited state. A layer of photoresist covers 
the substrate’s surface to capture the energy from the 
fluorescing excited state and the spin pattern can then 
be transfered to the photoresist. 

To see how MRL can achieve high resolution, let us first 
consider how to write a narrow line in a spin pattern. If 
we apply a gradient magnetic field, the Zeeman shift of 
the m = +/ — \ sublevels is a function of position. A 
narrowband microwave tt pulse is applied to excite the 
spins from m = 0 to one of the m = -|-/ — 1 sublevels. 
Only spins with the proper Zeeman shift are resonant 
with the microwave and the others are not appreciably 
excited. This is similar to the “resonant z-slice” in MRI 
[ 2 ^. The absolute resolution limit of this technique is 
given by [1^ 


Ax = 


r 

dAzM/dx 


( 1 ) 


where T is the spin transition linewidth and dAzM/dx 
is the spatial gradient of Zeeman shift. The width of the 
line is not limited by the wavelength but the spin transi¬ 
tion linewidth and the gradient of the magnetic field. In 
principle, we can achieve arbitrary resolution by increas¬ 
ing the gradient of the magnetic field. For example, if 
the dephasing time is about 10/rs [ 2 ^ and the gradient 
of the magnetic field is IQQGauss /the resolution can 
be a few nanometers. 

We can also address many spatial locations in parallel 
as in conventional MRI where the narrowband microwave 
is replaced by a pulse with more complex spectral con¬ 
tents . For example, if we want to print an arbitrary 
ID pattern characterized by /(x) whose value is normal¬ 
ized into 0 to 1 range, we should have excitation proba¬ 
bility given by [s^ 

fix) = ^{1 - cos[I2fl(x)Tp]} (2) 

where Tp is the pulse duration and flfl(x) = ^ab£{x)/h 
is the resonant part of the Rabi frequency seen by atoms 
at the position x. Here, £’(x) is the amplitude of the mi¬ 
crowave frequency component which is resonant to the 
transition frequency w(x) at position x, /iab is the transi¬ 
tion dipole moment between the two relevant energy lev¬ 
els, and h is the Plank constant. The transition frequency 
at position x is given by a;(x) = wq -b g^BB(x)/h where 
Wo is the zero field splitting, ggs = 2.8MHz/Gauss, and 
B{x) is the magnetic field at position x. The pulse in the 
time domain is the Fourier transformation of the required 
spectrum which in turn is related to position because of 
the magnetic gradient. 

2D patterns can also be printed in this setup where 
all three energy sublevels are involved and both X and 
Y magnetic gradients should be applied. The proposed 
pulse sequence is shown in Fig. 1(c). For example, we 
want to write a letter “T” (Fig. 1(d)). First, all the 
spins are polarized to the state | — 1) (green down trian¬ 
gles). This is the initial state. To write column 1, the 
X-gradient magnetic field is turned on and a narrowband 
microwave tt pulse (A), whose frequency is tuned to be 
resonant with spins in column 1, rotates spins in this 
column from | — 1) to |0) state. Then the X-gradient is 
turned off and the Y-gradient is turned on and another 
microwave pulse(B) with the desired spectrum rotates 
the spin in the first row from |0) to | -I- 1) (blue up trian¬ 
gles) state. To write another column without crosstalk, 
we need to turn on the negative X-gradient magnetic field 
again (and turn off the Y-gradient) and then apply the 
negative of pulse A to rotate unwanted spins (red circle 
dots) in column 1 back to the state | — 1) (green down 
triangles). Finally, we again turn on the negative Y- 
gradient (with X-gradient off) and the negative of pulse 
B is applied to rotate the spins in state |-|-1) back to state 
|0). A ID pattern is then formed in the column 1 on the 
|0) state. After that we can move to the next column and 
repeat the above procedures. As seen in the figure, the 
desired 2D pattern “T” can then be encoded in the |0) 
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state. Finally, an optical readout of the spin population 
can excite the spins in state m = 0 into a fluorescing 
excited state. If the photoresist is within the near-held 
region of the spins, instead of fluorescing to free space, 
most of the optical energy of the fluorescing excited state 
is absorbed by the photoresist through Foster resonance 
energy transfer (FRET) [3l|, H^. 

The FRET efficiency is given by e = [1 -|- 
where Rq is the critical FRET radius and R is the dis¬ 
tance between the emitter and the receiver. To capture 
the high resolution features of the patteren, the photore¬ 
sist needs to be within the FRET radius from the para¬ 
magnetic system. This could be achieved in several ways 
depending on the nature of the paramagnetic system. 
For example, for paramagnetic defects in crystals, this 
could be achieved by shallow implantation or engineered 
growth techniques. A more practical approach, however, 
would be to coat the substrate surface with a resin that 
contains paramagnetic molecules or nanoparticles with 
paramagnetic defects. 


III. NUMERICAL SIMULATION 


(a) 



t(|js) 



In the following, we apply a numerical simulation to 
write a ID and a 2D pattern onto the |0) state using 
MRL. The dynamics of the system is govern by the mas¬ 
ter equation 

'I 

p{x,y,t) = --[H{x,y,t),p{x,y,t)] - Cp (3) 

where p is 3 x 3 density matrix for the three energy sub- 
levels, H is the Hamiltonian which is given by 

Ai(a;) Di(t) 0 
Di(t) 0 D2(t) , 

0 Q.2{t) A2(y) _ 

where Hi (f) is the Rabi frequency of the | — 1) —>■ |0) tran¬ 
sition and Ai(x) is the corresponding detuning. D 2 (t) 
is the Rabi frequency of the |0) —> | + 1) transition 
and A 2 {y) is the corresponding detuning. The last term 
Cp = (po.iti + P±i.o)/’r 2 is the dephasing of the system. 
The whole process needs to be done before the population 
decay time Ti which can be of the order of milliseconds. 
Here we will neglect population decay but keep the de¬ 
phasing process. In the numerical simulation, we apply 
the Fourth-Order Runge-Kutta method to solve the time 
evolution of Eq. (3) with Rabi frequencies and detunings 
vary as the pulse sequence shown in Fig. 1(c). 

For the ID pattern example, we write a chirped grat¬ 
ing within a 150nm field of view using numerical simu¬ 
lation. Here, only an X-gradient and the | — 1) —|0) 
transition are involved. The spins are initially prepared 
in I — 1) state. We assume the wire diameter is l/rm and 
the distance between the centers of the two wires is 2pm. 
The currents in both wires are 50mA and in the same 
direction. The currents generate a magnetic gradient of 


H{x,y,t) 

h 



FIG. 2: (Color online) (a) Complex RF pulse used to write a 
chirped grating. Red line with cross symbols is the real part 
while blue line without symbols is the imaginary part, (b) 
The spatial distribution of the probability that the spins are 
prepared in the bright state, (c) The FRET efficiency on the 
substrate surface when the readout contrast is 100%. 


about 200Gauss/pm in the X-direction. The RF mag¬ 
netic field envelope of the microwave excitation pulse for 
the I — 1) —^ |0) transition in the time domain is shown in 
Fig. 2(a). The total time duration is about 20ps. This 
desired pulse could be generated with a commerical ar¬ 
bitrary waveform generator (AWG) [1^. Using an AWG 
also allows any non-uniformities in magnetic bias or gra¬ 
dient fields to be compensated by applying the appropri¬ 
ate inverse non-linearity to the AWG instruction set. In 
the numerical simulation, we assume the Rabi frequency 
has 10% amplitude fluctuation. The paramagnetic ma¬ 
terial is assumed to be 5nm deep. The probability that 
the spins are in the m = 0 state is shown in Fig. 2(b) 
where we can see that the desired pattern is encoded in 
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FIG. 3: (Color online) Numerical simulation to write the 
Texas A&M University logo into a two-dimensional area us¬ 
ing MRL. (a) The spatial distribution of the probability that 
the spins are prepared in the bright state, (b) The FRET 
efficiency on the substrate surface when the readout contrast 
is 100%. 


the TO = 0 bright state. The pattern of the internal state 
can then be transferee! to the photoresist by the FRET. 
The FRET efficiency on the substrate surface is shown 
in Fig. 2(c) where Rq is assumed to be 5nm. We can see 
that a chirped grating can be generated with very high 
contrast. 

For a 2D pattern numerical example, we write 
the Texas A&M University logo (“ATM”) into a 
72nmxl80nm region. This area is divided into 24 slices 
in the x direction and 60 slices in the y direction. The 
separation of the spins is chosen to be 3nm on a square 
grid. Each line consists three spins. The pulse sequence 
is shown in Fig. 1(d) where both X-gradient and Y- 
gradient magnetic fields are involved. In the numerical 
simulation, the gradients in both directions are generated 
by 100mA currents using the same wire separation as in 
the ID example. The narrow band microwave tt pulse 
used to drive transitions | — 1) ^ |0) is assumed to have 
a Gaussian time envelope. The time duration of the en¬ 
tire pulse A which selects a line is about 5.6/xs where the 
Gaussian envelope portion has a width of about 1.5/xs. 
The time duration of the complex pulse B is also about 
5.6^s which is discretized into 7000 time steps. The total 
time duration to write one line is about 22.4^s and it is 
about 500/rs for the whole process. The dephasing time 
is assumed to be about 20/iS [s^. Again, the amplitudes 


of the pulses have a fluctuation of about 10%. The prob¬ 
ability that the spins are in the bright state (to = 0) is 
shown in Fig. 3(a) where we can see that the desired pat¬ 
tern is successfully encoded into the internal state of the 
system. Finally, an optical readout can transfer the pat¬ 
tern of the internal states to the photoresist. Assuming 
that the paramagnetic systems are about 5nm deep and 
the FRET radius Rq is also 5nm, the relative FRET effi¬ 
ciency at the surface is shown in Fig. 3(b) where we see 
that the desired pattern can be formed. The full width 
at half maximum of the line in the middle is about 7nm. 


IV. POSSIBLE MATERIALS FOR MRL 

As is mentioned in Sec. II, MRL requires a solid sys¬ 
tem with high spin transition quality factors and high 
contrast optical readout of the spin state. There are 
a number of possible systems which can satisfy these 
requirements. One promising candidate is the negativ- 
ley charged Nitrogen Vacancy (NV) color center in di¬ 
amond. The NV exhibits an optically detected mag¬ 
netic resonance (ODMR) spin system with long coher¬ 
ence times, even at room temperature Coher¬ 

ence times larger than a millisecond have been experi¬ 
mentally achieved for NVs leading to their frequent use 
as nanoscale magnetic field sensors [s^. Such properties 
could give MRL resolutions as high as a few nanometers. 
Unfortunately at room temperature the optical readout 
of the NV spin state has less than 30% contrast. Nonethe¬ 
less the inherent nonlinearity of resist exposure may al¬ 
low us to produce high contrast patterns even with this 
relatively low optical contrast. 

One alternative candidate is the silicon vacancy (SiV) 
center in diamond. The SiV has an optical readout con¬ 
trast near 100% at low temperature since its ground 
state sublevel transitions can be optically resolved [s^. 
Unlike NV and many other materials, this high resolu¬ 
tion persists even for ensembles, due to the low inho¬ 
mogeneous broadening and low spectral diffusion of SiV. 
Granted that working at low temperature might produce 
some technical difficulties, nonetheless photoresist expo¬ 
sure has been demonstrated in cryogenic temperatures 
[40l| . Another important advantage of the SiV is that 
it can have optical Raman transitions connecting the 
ground state sublevels. Here light shift gradients could 
also replace magnetic field gradients [4l| so that no wires 
would be needed on the substrate at all. 

Finally we briefly consider other possible materials sys¬ 
tems for performing MRL at room temperature. One 
of these is the STl color center exhibiting a high op¬ 
tical readout contrast up to 50% albeit with shorter 
spin lifetimes [l^. Ruby is also of interest as it has 
shown room temperature optical Raman transitions with 
readout contrast comparable to NV[4^. Other alter¬ 
natives are the rare-earth doped materials. For in¬ 
stance, super-resolution microscopy has been demon¬ 
strated using praseodymium-doped yttrium aluminum 
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garnet (PriYAG) at room temperature [ij. The upcon- 
version UV emission of (Pr:YAG) is optimal for photore¬ 
sist exposure. Hence the prospects are good that a mate¬ 
rial suitable for room temperature MRL will eventually 
be identified. 


V. SUMMARY 

In summary, we have proposed an approach for optical 
lithography which is the inverse of magnetic resonance 
imaging. The technique uses atomic coherence in an en¬ 
semble of spin systems whose final state population can 
be optically detected. We design a pulse sequence which 
can write arbitrary ID and 2D patterns into the spin 
states. We presented a numerical simulation to prove 


the fundamental physics of MRL which shows that our 
method here is capable of producing arbitrary ID and 2D 
high-resolution patterns in one of the spin sublevels with 
high contrast. The spin pattern can then be transfered to 
the photorsist pattern by FRET. We also discuss several 
possible materials for our MRL method. 
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